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We describe theoretically a setup in which a tapered optical nanofibre is introduced into an 
optical lattice potential for cold atoms. Firstly, we consider the disturbance to the geometry of the 
lattice potential due to scattering of the lattice lasers from the dielectric fibre surface and show 
that the resulting distortion to the lattice can be minimized by placing the fibre at an appropriate 
position in the lattice. We then calculate the modifications of the local potentials that are achievable 
by transmitting off-resonant light through the fibre. The availability of such a technique holds 
the potential to deterministically create and address small well-defined samples of atoms in the 
evanescent field of the tapered nanofibre. 
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I. INTRODUCTION 

During the last two decades, advances in the eooling 
and trapping of atoms and ions has assisted in the cre- 
ation of clean and highly controllable systems, in which 
fundamental quantum mechanical experiments can be 
carried out with very low levels of noise. This has led 
to several breakthrough successes in the quest for imple- 
menting ideas of quantum information processing (QIP) 
[l[ , high-precision atomic clocks and quantum metrol- 
ogy Hi- 

For neutral atoms optical lattices have been important 
and hold a great deal of promise in this area. The high 
degree of control one has over the laser parameters has 
allowed for the execution of many seminal experiments 
in these periodic systems. In particular, by controlling 
the amplitude of the lasers one can adjust the trapping 
depth, which can act as a switch between regimes in 
which the dynamics are controlled either by tunnelling 
between different lattice sites or by interactions between 
the atoms. This has led to the celebrated observation of 
the superfluid-Mott insulator transition, in which a state 
with one atom per lattice site can be created 0, Q . 

States which have a well-defined number of particles, 
so-called atomic Fock states, are currently of large in- 
terest in physics. Their sub-Poissonian number statistics 
is valuable for applications in atom-metrology, quantum 
information processing and has merit for investigating 
the foundation of quantum mechanics as well. Several 
ground-breaking experiments have recently reported the 
creation of such states 04^ and a significant amount of 
theoretical work has been devoted to their characteriza- 
tion [lol.[T]|. Knowing the exact number of atoms a priori 
in each run of the experiment is still a difficult task, and 
techniques which can deterministically create a desired 
atom number are under vigorous development. 

Here we present a near-field optics approach to cre- 
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ating such definite atom-number states and propose the 
use of the evanescent field of an optical fibre as a tool 
for manipulating the optical lattice potential locally (see 
Fig. [1]). While standard optical fibres have a diameter 
of several hundred /Ltm, recent progress in tapering tech- 
niques allows for the creation of fibres of subwavelength 
diameter [l2[ and even down to 50 nm[l^. A significant 
amount of the intensity in these fibres is carried in their 
evanescent field and can therefore be used to create an 
optical potential for ultracold atoms. In this work we will 
examine the effects of introducing a sub-micrometer fibre 
into an optical lattice and demonstrate the possibility of 
deterministically creating states of fixed particle number 
using appropriately chosen fields inside the fibre. 

The paper is organized as follows: in Sec. |TT] we will 
present a short review of the potential forces relating 
to optical lattices and sub-micrometer diameter, single- 
mode silica fibres. We then discuss the modifications of 
an optical lattice potential due to effects from light scat- 
tering on a fibre in Sec. IIIII and examine several achiev- 
able potential geometries resulting from the combination 
of the lattice potential and the evanescent field poten- 
tials in Sec. HVl The resulting atomic state is discussed 
in Sec. fVland we finally conclude in Sec. 



II. POTENTIAL FORCES 
A. Optical Lattices 

To understand the influence the introduction of the 
fibre into an optical lattice has, let us first briefly re- 
view the optical potentials associated with optical lattices 
and nano-fibres. Optical lattices exist by today in many 
laboratories and represent periodic arrays of micro-traps 
generated by the dipole force of a standing wave laser 
light field [J H, • A variety of trapping geometries are 
achievable, with the most common being rectangular Q 
or triangular [Tsj . 

The simplest case of an optical lattice trapping poten- 
tial is given by a one-dimensional model, in which two 



counter-propagating laser beams interfere. This results 
in a standing wave for the optical intensity given by 



/(z) = /osin2(fcz) 



(1) 



where fc = 27r/A is the free space wave number of the laser 
light, Iq is the maximum intensity of the laser beam and 
the periodicity is given by A/2. The spatially varying ac 
Stark shift then forms a potential for the induced dipole 
moment, p, of the atom given by 

C/dip = -i(p-E) = -^K(a)/ (2) 

where eq is the vacuum permittivity, c the speed of light 
and a{uJL) the optical polarizability, which depends on 
the frequency of the laser field, E. Its real part is given 
by M 
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where the gj and ga are the statistical weights of the 
excited and ground states, respectively, the ujja are the 
transition frequencies and the jja are the emission tran- 
sition probabilities Depending on the detuning of 
the laser beam, the atoms can be forced to gather at the 
nodes or anti- nodes of the laser intensity pattern by using 
light blue-detuned (ujl > wq) or red-detuned {cul < wq) 
with respect to the chosen transition w, respectively. 

By introducing pairs of counterpropagating lasers in 
the remaining directions of space, higher dimensional lat- 
tices can be created. The interference terms between the 
laser fields in the different directions can be eliminated 
by choosing perpendicular polarization vectors of the two 
laser fields, which for a two-dimensional setup results in 
an intensity pattern represented by the sums of purely 
sinusoidal orthogonal fields (see Fig. [T|) 



= /q [sin^(fcx) 
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Throughout this paper we will consider this type of two- 
dimensional optical lattice, however a generalization to 
three dimensional, layered lattices is straightforward. We 
will also assume that every beam is independent and not 
created through retroreflection. 

Optical Lattices typically have lattice constants in the 
range of 400 nm to 650 nm and we will consider a lattice 
with a trapping wavelength A/2 = 527 nm. We choose 
the lattice to be loaded with ^'^•^Cs atoms, which localize 
in the high field regions. 



B. Sub-wavelength diameter optical fibres 

Recent developments in tapered, dielectric fibre tech- 
nology have made it possible to produce fibres with radii, 
a, as low as a few hundred nanometres (Tsj . In such fi- 
bres the core has vanished and they can be described 





FIG. 1: (Color online) Schematic showing a typical optical po- 
tential geometry resulting from the combination of an optical 
lattice potential and a trapping potential around an optical 
nano-fibre. 



by one large refractive index step between the remaining 
cladding, ni(a;), and the outside vacuum, n2. An inter- 
esting consequence of the subwavelength nature of the 
diameter is that the majority of the field will be guided 
in the evanescent field on the fibre's surface. It there- 
fore becomes accessible to atoms in the fibre's vicinity 
and light blue-detuned with respect to the atoms tran- 
sition frequency will create a repulsive force preventing 
the atoms from coming too close to the fibre which is at 
room-temperature. At the same time red-detuned light 
will result in an attractive force and a combination of 
both fields was suggested by Le Kien et al. [l^ as a way 
of creating a trapping potential around the fibre. This 
was experimentally observed in [Tsj . 

Let us briefly review the description of such a potential, 
following closely [l^. We consider two frequencies, 
and Wb, where the indices correspond to the red- and 
blue-detuned fields, respectively. They are chosen such 
that the single mode condition 



V, = h 
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is fulfilled [20[ and both light fields are in the funda- 
mental mode HEii. The intensity distribution of the 
evanescent fields depends on the polarization of the input 
fields and here we choose circular polarization for both 
beams to achieve angular symmetry [2ll . [2^ . In cylin- 
drical co-ordinates {r,(j),9}, the time-averaged intensity 
outside the fibre is then given by 

m^^e^ [K^ {q,r) + w,Kf {q,r) + f.Ki {q,r)] . (6) 

Here the Kn are the modified Bessel functions of the 
second kind and is the strength of the electric field. 
The decay of the fields from the surface of the fibre is 
characterized by g^, which is the reciprocal of the decay 
length Ai and given by 
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where i s the longitudinal propagation constant of the 
mode jm . Finally, the pre- factors are given by [22| 



Wi 



■2q^ 



(1-5.)^ ' 



(8) 
(9) 



3 



with s defined as 



III. SCATTERING AT THE FIBRE 
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and hi = {nfkf — (if)^. The combined optical potential 
around the fibre for a blue- and a red-detuned field is 
therefore given by (see Fig. [T]) 
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[Kl{qi,r) + WbKl{qi,r) + fbKHqbr)] 

[K^iqrr) + WrKliqrr) + /,.X|(g,r)] , 

(11) 

where the factors in front of the mode-structure terms 
are directly proportional to the powers of the individual 
light fields, Pr and Pb- 



C. Van der Waals interaction 

Finally we need to take into account the van der Waals 
attraction between the atoms and the fibre. The classical 
van der Waals potential felt by an atom near the surface 
of a dielectric fibre of infinite length is given by [1^ 

V{r)^-^ E / dk[eK^{kr) + {k^ + -)Kl{kr] 
47r'^eo ±^ Jo 
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Here the /„(a;) and Kn{x) are the modified Bcsscl func- 
tion of the first and second kind, respectively. It should 
be noted that this approximation neglects the resonant 
frequencies of silica. However, as these are substantially 
different and weaker than those of Cs atoms, this is jus- 
tified 0. 

A detailed analysis of expression ([T2|) was carried out 
by Le Kien et al. [l^, who found that for atoms close 
to the surface the van der Waals potential tends to the 
same values as that for a flat surface. The latter has the 
simple and well known form 
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The ground state Caesium atom has its dominant (-D2) 
line at 852 nm, which gives a van der Waals constant of 
C3 = 27r • 1.56kHz [23] and this value will be used 
throughout this paper. In the following we will also use 
the simplified expression (jl4[) whenever justified while 
making sure the full expression gives identical results. 



When a fibre is introduced in a position perpendicu- 
lar to the transverse plane of a two-dimensional optical 
lattice, the four incident beams will be scattered from 
the cylindrical surface and distort the regularity of the 
lattice. To describe this we approximate the fibre by an 
infinite cylinder of radius a = 150 nm oriented orthog- 
onally to the lattice vectors and assume that the waves 
undergo a linear scattering process [25j . In cylindrical 
co-ordinates the four incident waves of the conventional 
lattice have the form of plane waves 



Ei{0j) = Eo cxp [ikr cos{e - Oj)] , 



(16) 



coming from the angles 9j ~ 0, 7r/2, tt, 37r/2 and the total 
lattice field is that given by 
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Ei = W(E,(0)+E/(7r))^ 



(17) 

Assuming that the beams along the 7r/2-direction are po- 
larized parallel to the cylinder axis and the ones along the 
TT-direction orthogonal, the scattered field can be written 
as in [2^, where the respective polarizations are con- 
tained in Eq. 
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where the Hn' are Hankel functions of the first kind and 
the scattering coefficients are given by 

Jn{a)J^{ma) — mJn{ma)J^{a) 

(2) f2V 

Hn '{a)Jn{ma) - ■mJn{ma)Hn (a) 
^ _ mJn{a)J'n{ma) - Jn{ma)J'^{a) ^^^^ 
mHi^'^ {a)J'^ {ma) - Jn{ma)H^?' {a) ' 

From this the complete field follows as 



(20) 



Etot = J {Ej + eI) +{E^+Ei 



(22) 



Here we note that all of the square plots in this paper 
span an area of 4.2 jim x 4.2 /im. This 4.2 /im corre- 
sponds to approximately four optical lattice wavelengths 
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and thus an eight by eight grid of traps. The optical in- 
tensity in the vicinity of a fibre of radius a — 150 nm for 
a lattice constant of A/2 — hTJ nm is shown in Figs. [2] 
(a) and (c). One can see that if the fibre is located at 
a minimum of the optical intensity (Fig. [2ja)), the lat- 
tice structure is almost unaffected. Positioning the fi- 
bre at an optical intensity maximum on the other hand 
(Fig. HJc)) leads to noticeable disturbances, which, nev- 
ertheless, leave the basic lattice structure intact. Clearly 
larger fibres will lead to more scattering, however the 
numbers chosen here are well in reach of experimental 
possibilities. 

In general the scattered radiation propagates as a 
cylindrical wave and its intensity falls off as the inverse 
power of the radial distance. Since the energy flow is only 
in the planes of constant z, the scattered radiation cor- 
responding to a particular incident ray will be observed 
only in that plane which contains the incident ray and no 
scattering into other layers of a three-dimensional lattice 
occurs. 

The overall potential seen by the atoms must include 
the van der Waals potential and Figs. (Hb) and (d) show 
that the lattice sites most affected by the scattering are 
also strongly affected by the van der Waals potential 
(Note that the Cs atoms we are considering here are high 
field seekers). It is clear that in a shallower lattice the 
effect of the van der Waals attraction will be more severe 
on a larger range and we will show in the next section 
that the introduction of a repulsive blue field can be a 
useful tool for counteracting this effect. 
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FIG. 2: (Color online) (a) Lattice intensity including the scat- 
tering of the light on the fibre when the fibre is placed at an 
intensity minimum of the lattice, (b) van der Waals poten- 
tial and optical lattice potential which includes the scattered 
lattice field, (c) and (d) show the same for a fibre placed at 
a lattice intensity maximum. Each plot spans an area of 4.2 
^xm X 4.2 \xm and the lattice depth is 60 E_h. A colour scale 
has been used which varies from blue to red where the blue 
areas correspond to minima and red areas to maxima. 



IV. ADDING FIBRE POTENTIALS 
A. Compensating the van der Waals potential 

In order to minimize the disturbance of the lattice due 
to the van der Waals potential, we will study the possibil- 
ity of compensating the attractive potential with a repul- 
sive one from a blue-detuned optical field. The joint po- 
tential is simply given by adding the blue part of cq. (|lip 
to the van der Waals expression of cq. (|14p 




Since the modified Bessel functions have an exponentially 
decaying form it is not possible to perfectly compensate 
the van der Waals potential at all distances from the 
fibre. However, the discrepancy is weaker at larger dis- 
tances, which allows the reduction of the radius in which 
the attractive potential is significant. In Figs. [3] and |3] we 
show the potential for the two different positions of the 
fibre for different intensities of the blue beam. One can 
clearly see that in both cases it is possible to achieve a 
situation in which almost all lattice sites close to the fi- 
bre are still intact (central panel). This is important for 



the construction of a well-defined Mott insulator state 
around the fibre, which is a prerequisite to loading the fi- 
bre potential with a well-defined particle number. While 
the exact number of restored lattice sites also depends 
on the lattice depth, the graphs show typical achievable 
experimental values. For very small distances from the 
fibre surface, the attractive van der Waals potential will 
always be stronger than the compensating optical field 
and tunnelling will become an important loss factor at 
longer times. 

Let us also remark that with such a localized potential 
it is possible to remove atoms from specific lattice sites 
by using the fibre as a dark absorber [T^. Used in con- 
junction with existing ideas involving optical conveyor 
belts [l^l this simple setup could be an effective method 
of removing entire rows or patterns of atoms. 



B. Loading the fibre potential 

In the following we will consider the situation where 
an attractive, red-detuned field is added to the fibre as 
well. This will allow for the creation of a circular po- 
tential minimum around the fibre, deep enough to trap 
ultracold atoms. Recent experiments have demonstrated 
this by stochastically trapping atoms from a surrounding 
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FIG. 3: (Color online) Combined potentials (van der Waals, 
blue-detuned evanescent field and optical lattice) for a fibre 
placed at a minimum of the lattice intensity. The wavelength 
of the evanescent field is A;, = 700 nm and its power is in- 
creased through Pb =0, 0.05 and 0.10 mW from left to right. 
The lattice depth is chosen to be — 60i?rcc. 




FIG. 4: (Color online) Same as Fig. [3] with the fibre placed 
at a maximum of intensity of the lattice. Each plot spans an 
area of 4.2 fim x 4.2 ^m. 



thermal or condensed cloud [l8|. Since in our situation 
the environment around the fibre is given by the well or- 
dered optical lattice, a controlled melting of the lattice 
by the evanescent field will transfer a controllable number 
of atoms from the individual lattice sites into the fibre 
potential. The resulting state is therefore highly number 
squeezed and can be used in applications in quantum in- 
formation or metrology 0, [2^ . 

We study this process by assuming a realistic experi- 
mental situation of a Mott insulator made from Cs atoms 
with a resonant transition at Aq = 852 nm Q. For the 
two light fields in the fibre we consider a blue-detuned 
field at a wavelength of Af, = 700 nm and the red- 
detuned at Xr — 980 nm. The dctunings of the fibre 
fields from the dominant line of the atom are then given 
by ^ = -46 THz and #^ = 76 THz and with a fibre 
radius of 150 nm, the evanescent decay lengths corre- 
sponding to the blue and red fields are Ab = 0.36 /xm 
and Ar = 1.83 /im. The two-dimensional optical lattice 
we consider has a depth QOEr. 

Two typical examples of resulting trapping geometries, 
when all potentials are taken into account, arc shown in 
Figs. [5] and [6] for a fibre located at an intensity minimum. 
Since we assume the fibre to be initially dark, the lattice 
sites which are visible closest to the fibre in Figs. [SJc) 
andini^c) are actually empty due to being within the ra- 
dius of the surfaces van der Waals potential. For all other 
sites in the vicinity of the fibre one can clearly see that 
the addition of the red and blue fields allows for the lower- 
ing of the on-site energies. Therefore, as the asymptotic 



(a) 



(b) 2 

1 

MK 

-1 

-2 







4.2 pm X 4.2 Mm 



\(c) 




2 2 
Distance in units of optical lattice 
wavelength 



FIG. 5: (Color online) (a) The black solid line represents the 
resulting two-dimensional potential geometry if the red and 
blue detuned beams in the fibre are switched on for Pb = 0.12 
mW and Pr — 0.036 mW. The lattice sites for which the 
energy is lowered such that the residing atoms will be trapped 
in the fibre potential after the optical lattice is switched off are 
marked with a cross, (b) and (c) cut through the potential 
at the lines indicated in (a). The relevant on-site energies 
are indicated as well and the dotted blue line indicates the 
potential after the optical lattice is switched off. The fact 
that the fibre potential does not always meet exactly the zero 
points of the overall potential before switch-off is due to phase 
shifting of the optical lattice after scattering on the fibre. 



potential of the evanescent field goes to zero and since 
there is no local maximum in the fibre potential, a sud- 
den switch off of the optical lattice will leave all atoms 
with centre-of-mass energies < trapped in the fibre po- 
tential alone. In the example shown in Fig. [SI where 
Pb = 0.12 mW and = 0.036 mW, one finds that this 
condition is fulfilled for 8 lattice sites. Increasing the 
evanescent fields to Pb = 0.17 mW and Pr = 0.082 mW 
(see Fig. [5]) the radius of the evanescent field increases 
and careful examination shows that 20 sites are reached. 
Note that, due to the rectangular geometry of the consid- 
ered lattice, only certain atom numbers can be achieved 
and realistic parameters limit this technique to samples 
of only a few tens. If the switch-off process of the optical 
lattice is done on a time-scale shorter than the typical 
atom tunnelling time in an optical lattice (which is of 
the order of several milliseconds @ ) , no other atoms will 
have been able to join the well-defined sample. 

The state created in this sudden-switch off procedure 
is a highly excited, out-of-equilibrium state in the final 
potential and subsequent cooling is necessary to prevent 
further atom loss due to scattering and thermalisation. 
While a detailed calculation of these efi'ects go beyond the 
scope of this work, loss through re-thermalisation can be 
minimised by applying a Fcshbach resonance while carry- 
ing out this process and subsequently adiabatically low- 
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FIG. 6: (Color online) Same as Fig. \5\ but for higher powers 
of the evanescent fields, Pt = 0.17 mW and Pr = 0.082 mW. 
A larger number of atoms will be trapped after the optical 
lattice is switched off. 



ering the fibre potential while switching the interaction 
back on (30| . 

It is also worth pointing out that in both examples 
above the presence of the repulsive blue field ensures the 
existence of a repulsive wall between the fibre and the 
atoms, thereby preventing direct atom loss through the 
room temperature object. However, the power in the 
blue detuned field does not correspond to the same field 
strength that optimally cancels the effects of the van der 
Waals potential, as discussed in the previous chapter. It 
is rather necessary to over-compensate the van der Waals 
potential and recreate the trapping minimum using the 
red-detuned field. In the next section we will discuss the 
nature of the atomic state created in the fibre potential. 



V. ATOMIC STATE 

Let us finally briefly characterise the atomic many- 
body state that can be created by the procedure above 
and focus in particular on the nature of the correlations in 
the sample. For this we first consider the effective dimen- 
sionality of the ground state of the potential around the 
fibre after the optical lattice is switched off. Since the size 
of the radial ground state of the potential will be much 
smaller than the curved, azimuthal one, we can assume 
an approximate separation of the wavefunction in the two 
directions. This allows us to describe the spectrum in 
the azimuthal direction by a free-space periodic potential 
with the well known spectrum = (n^?i^7r^)/(2mL^), 
where L = 27rrm is the circumference of the potential 
at the position of its radial minimum, ?■,„. Since no 
analytical expression for the position of this minimum 
is known, we find it numerically and estimate the en- 



ergy difference between the ground and the first excited 
state in the azimuthal direction to be of the order of 
AE^Q = E^ ~ E^ ^ 10^3^ J, for both situations shown 
in Fig. [S] 

To find the spectrum in the radial direction, we nu- 
merically diagonalize the radial part of the potential for 
a wide range of parameters and find typical values for 
the separation of the ground and first excited state to be 
of the order of A_E[q ^ 10~^° J. This significant differ- 
ence in the stiffness of the spectra in the two different 
directions (the z-direction can be adjusted separately to 
be stiff) translates into an approximate one-dimensional 
situation with an aspect ratio of ^ 10'^. 

Having established the effective dimensionality of the 
potential, the many body state of a one-dimensional 
Bose gas can now be characterized using the Lieb- 
Liniger parameter 7 = mgin /fi^niz) [3l|. Here qid is 
the one-dimensional coupling constant given by gio = 
^^^{ai_ - Caac)"^ with C « 1.4603 and nio is the 
linear density of the atoms [13 . For values of 7 ^ 1 the 
atomic many body state would be in the strongly cor- 
related Tonks-Girardeau regime, whereas for 7 < 1 the 
gas can be treated as weakly correlated. For both cases 
shown in Fig. [5] the radial ground state size is of the order 
of aj^ ~ 0.2 /im (we assume the same is achieved in the 
z-direction) and the position of the radial minimum is at 
rm ^ 6.6 yum. This leads to values of 7 = 0.560 for the 
8 atom case and 7 = 0.346 for the 20 atom case, putting 
both states firmly in the weakly correlated regime. 



VI. CONCLUSION 

In this work wc have suggested that the combination 
of optical lattices and tapered optical nano-fibres can be 
used to create small atomic samples which allow control 
over the final atom number. While introducing the fibre 
into the optical lattice inevitably leads to a disturbance 
of the lattice in the vicinity of the fibre due to scatter- 
ing of the lattice beams, we have shown that this can be 
minimized and, due to the small fibre diameter, usually 
leaves the overall lattice structure intact. The attractive 
van der Waals potential close to the surface of the fibre 
can be compensated by using a blue-detuned evanescent 
field around the fibre, which allows a reduction of the 
range of the fibre's influence to the size of a single lattice 
site for typical experimental parameters. Adding a sec- 
ond, red-detuned light field to the fibre then allows local 
melting of the optical lattice and can be used to create a 
small sample with a well-defined atom number. Finally, 
we have shown that these samples are in the superfluid 
regime and therefore good candidates for investigating 
the physics of persistent currents, or, using more varied 
optical potentials around the fibre, the physics of super- 
fluid squids. 
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